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Analysis of oscillations of atmospheric neutrinos 
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We briefly review the current status of standard oscillations of atmospheric neutrinos in schemes with two, 
three, and four flavor mixing. It is shown that, although the pure v^^ -^ Vr channel provides an excellent 2v flt to 
the data, one cannot exclude, at present, the occurrence of additional subleading fp ^ Ve oscillations (J,v schemes) 
or of sizable Vp. — > Vs oscillations {J^v schemes). It is also shown that the wide dynamical range of energy and 
pathlength probed by the Super-Kamiokande experiment puts severe constraints on nonstandard explanations of 
the atmospheric neutrino data, with a few notable exceptions. 
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1. Introduction 

It is well known that the Super-Kamiokande 
(SK) atmospheric neutrino data can be beauti- 
fully explained in terms oi 2i> oscillations in the 
j/^ — > i/^ channel |]l|. This interpretation is also 
supported by MACRO 
Conversely, pure i^^ -^ i 
vide a good fit to the SK data Q, and are inde- 
pendently excluded by the negative i^e disappear- 
ance searches in the CHOOZ [|| and Palo Verde 
|6| reactors. Pure ly^ -^ Vs oscillations {us being a 
hypothetical sterile neutrino) are also disfavored 
by SK |i|,§ (and by MACRO §), due to nonob- 
servation of the associated matter effects [p[ and 
neutral current event depletion |b|. 



Although two-flavor v^ 



oscillations repre- 



sent the most economical explanation, it should 
be stressed that, to some extent, additional os- 
cillation channels may be open, as naturally ex- 
pected in iv and Av schemes pO| | accommodating 
the current phenomenology. Moreover, i'^ disap- 
pearance might be driven by dynamics different 
from the simple mass-mixing mechanism. In this 
article, we briefly review the status of such so- 
lutions, with emphasis on: (i) scenarios involving 
more than two states {iv and Av mixing) , and (ii) 
scenarios involving nonstandard dynamics (decay, 
extra dimensions, decoherence). 



2. iv oscillations 

Assuming that two out of three active i^'s are 
almost degenerate (say, mi 2± TO2), it can be 
shown [y that atmospheric z/'s probe only m? = 
m§ — m\ 2 and the mixing matrix elements f/as: 

3 i/ parameter space = (m^, C/^g, C/^3, t/^3) , (1) 

with J7g3 -I- C/?^3 + C/^3 = 1 for unitarity. 



Progress in m^ bounds for unconstrained 3i/ mixing 
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Figure 1. Progress in bounds on m^, as derived 
by iv analyses of atmospheric and reactor data, 
both before and after SK and CHOOZ. 
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Figure 2. SK zenith distributions, normalized to no-oscillation expectations. Dots with error bars: SK 
data. Dashed and sohd lines: best fits to SK only and to SK+CHOOZ pit]. 



We present a preliminary update lO] of previ- 
ous limits [Q on such parameters, using the lat- 
est data from SK (70.5 kTy) and CHOOZ |]. 
The SK data include 55 zenith bins: lO-flO bins 
for the subGeV (SG) e+fji events, 10+10 bins for 
the multiGeV (MG) e+/i events, and 5+10 bins 
for the upward stopping (US) and through-going 
(UT) /i events. For CHOOZ, we use the total rate 
(one datum). We accurately calculate all such ob- 
servables, and x^-fit them (see H] for details). 

Figure 1 shows that the allowed range for m^ is 
stable around 3 x 10"'^ eV^. The same figure also 
shows the impact of SK and CHOOZ in sharpen- 
ing [|ll],|| prior bounds on m? |l2| . 

Figure 2 shows the SK data and the best-fit 
theoretical distributions. The best fit for SK data 
only (x^ — 47.5, dashed line) is found at 



{m', [74, C/^^3, C/4) ^ (3.5, 0.07, 0.57, 0.36) , (2) 

where [m^] = lO'^ eV^. For t/f^ = 0.07, the 
theoretical MGe distribution shows a distortion 
which, however, is well within the uncertainties. 
The weak preference for U^^ 7^ is suppressed by 
CHOOZ data. The SK+CHOOZ best fit (x^ = 
49, solid lines) basically corresponds to pure i^^ -^ 
Vr oscillations with maximal mixing, 



with limited allowance for extra v^ mixing |ri| , 
SK data only : UJ^ < 0.31 (0.38) , (4) 



SK + CHOOZ : U^^ < 0.04 (0.07) 



(5) 



the bounds being at 90 (99%) C.L. for 3 d.o.f. 
Unfortunately, it appears very difhcult to probe 
(through present atmospheric data) values of C/gg 
as small as a few %, which may entail interesting 
Earth matter effects [MESJ. Constraining C/^g is 
a major task for future atmospheric |14|, reactor 
p^ and accelerator [|l6| v experiments. 

The bounds on Sv mixing are more evident in 
the {ue^Vfj^^Ur) triangle plot, embedding the uni- 



{m 






C/^g) ^ (3.0, 0, 0.5, 0.5) 



(3) 



tarity constraint (see p 12| for details). Figure 3 
shows the allowed regions in such triangle, whose 
lower and right sides represent the subcases of 
pure i^^ -^ Vr (allowed) and pure v^ -^ i/g (ex- 
cluded). Large v^ mixing is allowed by SK alone, 
but not by the SK+CHOOZ combination, where 
only a narrow region survives near the lower side 
of the triangle. In such region, [/^g ~ U"^^ within 
a factor of two [e.g., (C/^^, V^^) ~ (2/3, 1/3) is 
also allowed]. In conclusion, the iv analysis of 
SK+CHOOZ shows that the v^ — > i/g channel 
might be open with a few % amplitude. Future 
atmospheric, reactor, and accelerator v experi- 
ments will test this interesting possibility. 
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Figure 3. Three-flavor analysis in the (fe, 1^11, Vt) triangle plot, for five representative values of rn^. Left 
and middle column: separate analyses of Super-Kamiokande 70.5 kTy data and CHOOZ final data, 
respectively. Right column: combined SK+CHOOZ allowed regions. The SK+CHOOZ solutions are 
close to pure Vn ^^ Vr oscillations, with upper limits on t/| 



^^3 in the few percent range |11[ 



3. 4;/ oscillations 

The current evidence for v oscillations coming 
from solar, atmospheric, and LSND data can be 
accommodated by introducing a fourth, sterile 
neutrino state v^ |10[ . The mass spectrum seems 
then to be favored in the "2+2" form (two sepa- 
rated doublets) ||I^, although the "3+1" option 
(triplet plus singlet) is not dismissed llq|. 

In 2+2 models, it is often assumed that atmo- 
spheric V oscillations involve either the i/^ -^ Vt 
or the v^ — > v^ channel. Correspondingly, it is 
assumed that solar v oscillations involve either 
the i^e — > Vs or the v^ — > Vt channel. Such sim- 
plifying assumptions are challenged by the most 
recent SK data [|l|Jl9f|, which disfavor oscillations 
into Vs for both atmospheric and solar neutrinos. 
However, it should be realized that atmospheric 
^'^'s and solar I'e's may also oscillate into linear 
combinations of i^s and i^t- ||20|| (rather than into 
Us and Ur separately), e.g.. 



Atmospheric i^^ — > i/+ oscillation fit 



atm. neutrino oscillations 
solar neutrino oscillations 



where 
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with ^ to be constrained by experiments. A re- 
cent analysis of Ue — > t^- solar oscillations shows 
that all the usual solutions (MSW or vacuum) are 
compatible with solar data for sin^ S, > 0.3 pit]. 

Concerning atmospheric i^'s, we have analyzed 
lH the same data as in Fig. 2 for ^ G [0,7r/2]. 
Figure 4 shows some representative results of the 
X^ fit, as a function of the mass square difference 
TO^. The fit for unconstrained i^ (thick solid line) 
is almost equal to the one for ^ = (pure u^ -^ 
Vt, thin solid line), implying that the SK data 
prefer small or zero admixture of Vg- The case 
^ — tt/2 (pure f^ -^ Vg, dashed line) leads to 
Ax^ ~ 15 and is disfavored. However, the case 
^ = 7r/4 (fifty-fifty admixture of Ut and v^ , dotted 
line) leads only to a modest increase in x^ and 
cannot be excluded. 

The Av analysis can also be done in a trian- 
gle plot (different from the Zv case) embedding 
the (i/^, z^s, j/t) unitarity constraint p2[. Figure 5 
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1/+ = {v^+v,)/\i2 
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m (X 10"'eV') 

Figure 4. x^ fit of 55 SK data bins (70.5 kTy) for 
z/p — !■ I/+ oscillations, under various assumptions 
for the Vg component oiv^. A large Vg component 
(e.g., 50%, dotted line) is not excluded [p2|. 



shows the results for separate and combined SK 
data sets. It can be seen that, in the combination, 
the case of pure i/^ — > Vj- oscillations (left side) is 
allowed, while the case of pure i/^ — > Vg oscilla- 
tions (right side) is significantly disfavored. How- 
ever, there are intermediate solutions for sin^ £_ < 
0.7 which have a significant admixture of Vg. Such 
results and constraints emerge from the interplay 
of low-energy data (which are more sensitive to 
m?) and high-energy data (more sensitive to the 
Vg component through matter effects, scaling as 
i neutron density x sin^ S, |2^ ) . 

A qualitative comparison between such results 
P4 and those in ||21] indicates that atmospheric 
1/ data can be reconciled with any of the oscilla- 
tion solutions to the solar i/ problem in the range 
0.3 < sin^ ^ < 0.7. A somewhat different Av anal- 
ysis ||23| derives similar conclusions. Summariz- 
ing, it turns out that world v oscillation data are 
consistent with Ai/ solutions to the solar and at- 
mospheric anomalies, involving oscillations into 
both active and sterile states at the same time. 
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Figure 5. A,v analysis in the (t/^, z^s, i^t) triangle plot, for five representative values of rr? . First three 
columns: separate analyses of SGe+SG//, MGe+MG/i, and US/i+UT^ data. Right column: all SK data 



(70.5 kTy). The allowed regions typically include pure Vy^ 



oscillations (left side of the triangle) and 



disfavor pure Vu 



oscillations (right side of the triangle). However, intermediate situations with Vy^ 



mixing with both v^ and Vg are allowed inside the triangle 



4. Nonstandard dynamics 

The SK data probe three decades in pathlength 
L and four decades in energy E. Such a wide dy- 
namical range severely constrains deviations from 
the standard L/E behavior of the P^r transition 
probability, which are expected in the presence of 
exotic dynamics p4] (e.g., violations of relativity 



principles |25 2q], which lead to a L-E behavior). 

An analysis of older data (45 kTy) has shown 
that, assuming a L • £"' dependence of the phase, 
the SK measurements constrain n to be very 
close to —1, thus favoring standard oscillations, 
and excluding several nonstandard explanations 
P7|] . Such results, shown in Fig. 6, have been 
strengthened by the latest SK data j|l|. A pe- 
culiar FCNC scenario with n = (28| is also 
strongly disfavored — as any energy independent 
mechanism for f^ disappearance — by combining 
low and high energy SK data |£9[. Therefore, F^r 
seems to be (dominantly) a function of L/E. 

However, is P^t necessarily a periodic function 
of L/El The answer is, surprisingly, no. There 
are (at least) three exotic scenarios which predict 
a monotonic decrease of the oscillation probabil- 
ity in the relevant L/E range, and that are any- 
way reasonably consistent with the data. 

Bounds on n for unconstrained a and B 



P„T = a sir' I (3 



10' km GeVy 
SK 45 kTy: SGe,;i -^ MGe,M + UP/^ (30 data) 




n - -0.9 ±0.4 

g: 90% C.L. (N„f = 3) 



-2 -1 

energy exponent n 



The first scenario involves ly decay ||30|] , with a 
decay length of the order of the Earth radius. The 
second scenario [ pl| assumes j^^ mixing with neu- 
trino states propagating in large extra dimensions 
P4 |. A third scenario Q assumes nonstandard 
Liouville dynamics |34| , leading to v decoherence 
and thus to a damping of oscillations. Figure 7 
shows that the best fit for pure decoherence does 
not differ significantly from the standard oscil- 
lation one [Q. The two cases shown in Fig. 7 
correspond to different functional forms for Pfj,fj., 



oscillation : P, 



fj.fl 



decohere. : P„ 



^[l + cos{+pL/E)], (9) 
^[1 + eM-pL/E)] , (10) 



with [E] = GeV, [L] = km, and p ~ 7 x IQ-^ 
GeV/km in both cases. Such forms have the same 
asymptotic behavior, namely, {PfLfi) — 1(2) ^^r 
small (large) L/E, but they significantly differ 
for intermediate values of L/E where, however, 
the large energy-angle smearing of SK prevents a 
clear discrimination. 

Although such nonstandard explanations | |30| , 
|3l| , |33[ of SK data do not survive Occam's ra- 
zor, they survive the current experimental tests 
for a simple reason: the oscillation pattern (ap- 
pearance of Vr and re-appearance of v^) has not 
been directly observed so far, and a monotonic 
j/^ disappearance is not excluded yet. Therefore, 
the unambigous observation of an oscillation cy- 
cle represents an important task for future atmo- 
spheric Q and accelerator [|6| v experiments. 
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Figure 6. Bounds on the energy exponent n, 
assuming oscillation phase oc L ■ E^. (Older 45 
kTy SK data used in this figure [EtI.) 



Figure 7. Comparison of standard oscillations 
(solid line) vs neutrino decoherence (dashed line) 
as explanation of the SK data. See pSl for details. 



5. Conclusions 



Rujula, M.B . Gavela, and P. Hernandez, hep- 
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simple and beautiful explanation of the SK data 
(as well as of MACRO and Soudan2). How- 
ever one cannot exclude, in addition, subleading 
j/^ -^ Vf. transitions (possible in "iv models) or siz- 
able v^^ -^ Vs transitions (possible in Av models). 
Moreover, the nonobservation of an oscillation cy- 
cle still leaves room for exotic dynamics. Further 
experimental and theoretical work is needed to 
firmly establish both the flavors and the dynam- 
ics involved in atmospheric v^ disappearance. 
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